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RADIATION  TRANSPORT  IN  A  NITROGEN  PLASMA 
GENERAL  FORMALISM  AND  1-DIMENSIONAL  MODEL 

I.  INTRODUCTION 

An  energetic  electron  beam  propagating  through  air 
dissociates  and  ionizes  N2  and  O2  leaving  behind  a  hot 
plasma  composed  primarily  of  N,  N+,  N++,  0,  0+,  and  0++. 
Emissions  from  this  plasma  cool  the  heated  channel  as  well 
as  provide  diagnostic  information.  To  properly  describe  the 
evolution  of  the  channel  under  nonequilibrium  conditions  one 
must  solve,  simultaneously,  the  time-dependent  rate 
equations  which  provide  the  population  densities  of  all 
bound  and  continuum  states  involved  in  the  radiative  and 
collisional  processes  and  the  corresponding  radiative 
transfer  equations  which  account  for  the  effect  of  the 
important  emissions  on  the  plasma.  In  this  paper  we  show 
the  effect  of  radiation  transport  on  the  cooling  of  a 
nitrogen  plasma.  The  nitrogen  plasma  consists  of  N,  N+,  and 
N++ ,  each  of  which  exists  in  various  excited  states  as  well 
as  ground  states.  The  radiation  transport  is  treated  as  1- 
dimensional . 

Depending  on  the  nature  of  the  specific  problem,  there 

are  many  approaches  to  solving  the  radiative  transfer 

equation1-.  For  example,  geometry,  time-dependence, 

radiation  characteristics  ( line  versus  continuum  radiation, 

single  lines  versus  many  lines,  etc.),  optical  thickness, 

inclusion  of  scattering,  and  the  thermodynamic  state  of  the 

gas  all  play  a  role  in  determining  the  best  approach  and 

corresponding  set  of  approximations.  For  instance,  rarely 

is  an  exact  solution  attempted  for  geometries  other  than 
2 

planar  .  In  many  cases  decomposition  of  the  radiative 
transfer  equation  into  a  set  of  moment  equations  with 
corresponding  restrictions  on  the  angular  dependence  of  the 
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intensity  has  proven  to  be  a  useful  procedure.  Kulander 
used  such  an  approach  to  investigate  radiative  transfer  in  a 
high-temperature  nonequilibrium  nitrogen  plasma  confined  to 
a  planar  geometry.  For  radiative  transfer  in  optically 
thick  regions,  i.e.  regions  which  are  large  compared  to  the 
optical  depth  of  the  radiation,  the  diffusion  approximation 
method,  as  developed  in  neutron  transport  theory^,  may  be 
applicable.  Comparisons  between  moment  method,  diffusion 
method,  and  exact  (or  nearly  exact)  results  have  been  made 
for  model  problems^.  Multi-frequency  grouping  methods  are 
often  applied  when  it  is  necessary  to  transport  many 
frequencies  and  possible  to  define  groups  of  frequencies 
within  which  the  absorption  coefficients  are  similiar. 

7 

Average  radiative  transfer  quantities  are  then  calculated  . 
Combinations  of  the  above  approaches  may  be  undertaken  as 

O 

well  .  One  method  that  has  seen  extended  use  is  the  escape 

probability  formalism,  as  based  on  the  approach  first 

9 

introduced  by  Holstein  .  This  method  has  been  used  for  a 
variety  of  problems^  '  ^  . 

The  approach  presented  in  this  paper  derives  from  the 
escape  probability  formalism.  Spectral  line  transfer  is 
emphasized.  The  basic  assumption  of  the  theory,  therefore, 
is  the  absence  of  correlation  between  the  frequencies  of 
emitted  and  absorbed  photons  in  the  sequential  process. 

This  is  rigorously  satisfied  when  the  line  broadening  is 
determined  by  impact  mechanisms  as  opposed  to  the  Doppler 
effect.  Rather  than  attempting  to  solve  the  transfer 
equation  f requency-by-f requency  over  a  line-width,  an 
effective  probability  is  derived  for  the  propagation  of  any 
photon  whose  origin  is  a  single  transition.  This  approach 
allows  for  the  treatment  of  many  lines  and  accounts  for  the 
phenomenon  that  photons  emitted  from  line  wings  have  longer 
absorption  lengths  than  those  emitted  near  line  center.  It 
also  provides  for  a  description  of  such  well-known  radiative 
transfer  effects  as  channel  cooling  due  to  radiative 
emission,  diffusion  of  radiation  energy,  and  photo-pumping, 
i.e.  photoabsorption  leading  to  an  enhanced  production  of 
electrons,  ions,  and  excited  atoms/ions.  The  method  is 
valid  for  nonequilibrium  conditions  and  when  expected 
reduces  to  equilibrium  results. 


The  details  of  the  above  approach  and  its  relationship 
to  solving  the  chemistry  rate  equations  are  discussed  in 
section  II.  Specifics  of  the  model,  in  particular,  line 
profile,  geometry,  plasma  composition,  and  temperature  and 
energy  equations  are  presented  in  section  III.  Numerical 
results  for  two  models  with  specific  initial  conditions 
given  and  discussed  in  section  IV.  Concluding  remarks  are 
reserved  for  section  V. 
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GENERAL  FORMALISM 


II . 


A. 


Population  Dynamics 
The  coupled  rate  equations  which  describe  the  evolution 
of  the  nitrogen  plasma  include  the  following  collisional  and 
radiative  processes:  electron  impact  excitation  and  de¬ 
excitation  of  N  and  N+,  electron  impact  ionization  of  N  and 
N+,  three-body  recombination  of  N+  and  N++,  radiative 


recombination  of  N+  and  N++,  spontaneous  emission, 


stimulated  emission,  and  absorption  for  N  and  N+ .  The  rate 


.2-1 , 


equation  for  Nffl  '1'(r,t),  i.e.  the  population  density  of  the 


(z-1)  ion  (or  neutral  if  z-1 )  in  the  electronic  state  whose 
index  is  m,  is 
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In  Eq.  (1)  all  atom,  ion,  and  electro  population 
densities  are  implicit  functions  of  r  and  t.  The  rate 
coefficients  for  excitation,  de-excitation,  ionization, 
three-body  recombination,  radiative  recombination,  and 
spontaneous  emission  are  designated  by  X,  Y,  S,  ,  and 

A,  respectively.  Subscripts  refer  to  the  states  involved  in 


the  particular  transition  while  superscripts  identify  the 

2  —  1 

ion  stage.  For  example,  Xlm  is  the  rate  coefficient  for 
collisional  excitation  of  ions  with  net  charge  z-1  to  an 


electronic  state  m  from  state  1.  The  effects  of  radiation 


transport  on  individual  species  are  contained  in  the  last 
four  terms  which  account  for  stimulated  emission  and 
photoabsorption.  These  terms  depend  on  the  normalized 


frequency  profiles  for  absorption  and  emission,  $  and  , 


the  appropriate  Einstein  B-coeff icient,  the  linewidths  Av, 
and  the  mean  angle-averaged  intensity,  J  (r,t).  The  mean 
intensity  Jv  is  obtained  by  solving  the  radiative  transfer 
equation  for  frequency  v;  this  is  discussed  in  detail  later 


below. 


The  radiation  terms  in  Eq .  (1)  are  simplified  by 
assuming  the  emission  and  absorption  profiles  are  identical 
Using  the  standard  relations  between  the  Einstein 
coefficients  Eq.  (1)  becomes 
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Stimulated  emission  is  now  treated  as  a  modification  to 


absorption . 


■vV/ 
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If,  in  Eq.  (2),  a  "gross  re-absorption  factor"  is 
defined  for  each  transition  line  as 
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then  the  rate  equation  becomes 
dN2-1 

—  (r,t)  -  collisional  terms  +  spontaneous  emission 

dt 
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Optically  thick  and  thin  results  are  frequently  obtained  by 

treating  the  set  of  {A}  as  constants  with  values  ranging 

between  0  and  1.  In  this  manner  the  radiative  transfer 

equation  need  not  be  solved  and  limiting  population 

12  —  16 

densities  are  obtained 


B.  Radiative  Transfer  Equation 

The  radiative  transfer  equation  in  its  simplest  form  is 


ds 


(5) 


Here,  is  the  absorption  coefficient  (corrected  for 
stimulated  emission)  at  frequency  v,  j  the  emission 
coefficient  for  the  radiation,  I  the  intensity,  and  s 

defines  a  ray  through  the  medium.  The  general  non  steady- 
state  solution  to  Eq .  (5)  is  written  in  integral  form  as 


I„(r,t,0,T) 


s 


ds'  j  (s' , t' )  e 


Js-  ds” 


(6) 


In  Eq.  (6)  t'  -  t-(s-s')/c,  t''  =*  t-(s-s'')/c,  and  j  ds' 
the  amount  of  radiation  born  at  s',t'  on  ray  ds'  in  a 
direction  0, Y.  Multiplying  by  the  exponential  gives  the 
fraction  of  radiation  that  arrives  at  s,t.  The  coordinate 
system  is  shown  in  Fig.  1.  In  this  study  continuum 
absorption  and  emission  are  neglected.  The  emission 
coefficient  for  a  specific  transition  line  is 


.  z  —  1  .  ..z  —  1 

A  ,  <f>  N 

ml  v  m 


The  absorption  coefficient  for  the  same  transition  is 


_  2  —  1 
*1 


where 


ym  .,2-1  .z-1  , 
z-1  1  ml  v 


The  statistical  weights  for  the  upper  and  lower  states  are 
z-1  z-1 

and  ’i  • 


It  is  the  mean  intensity  J  (r,t)  that  enters  into  Eq . 

(2)  and  determines  the  effect  of  radiation  transport  on  the 

individual  population  densities.  Letting  p  *  |r  -  r'|, 

2 

defining  the  volume  element  dV  =  dr'  =  o'  sin  0  d0  dT  dp', 
and  averaging  over  0,T  the  mean  intensity  is 


where  t'  -  t-p'/c  and  t ' '  =  t-p''/c.  Eq.  (10)  gives  the 
angle-averaged  solution  to  the  radiative  transfer  equation 
for  a  specific  frequency  \> . 


C.  Population  Dynamics  Plus  Radiative  Transfer 

The  rate  equations  may  be  rewritten  in  simpler  form  as 


dNm  1 

—  (r,t)  =  collisional  terms  +  spontaneous  emission 

dt 
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(11) 


where  the  set  {R}  is  defined  by  comparison  with  Eq .  (2). 
Substituting  Eq.  (10)  into  the  absorption  term,  noting  that 
is  given  by  Eq.  (8),  and  interchanging  the  \>  and  r' 
integrations  results  in 


Chr.t) 


dc'  Nfhc.f) 


’ml 
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(  12) 


Gm^  is  the  probabilty  that  a  photon  emitted  at  r',t'  is 
absorbed  in  a  volume  element  dr  at  r,t  and  is  given  as 


Gml(r''t''r't) 


4  ii 


r-r'l2 


dv  4>v  k^(  c,  t )  e 


o  dp" 


.  (13) 


ml 


The  discretized  version  of  Eqs.  (11)  -  (13)  forms  the  basis 
for  the  "escape-probability  formalism"  discussed  previously. 

It  is  useful  to  compare  Gm^  to  Tmj_  (  p  '  ;  t )  ,  i.e.  the 
probability  a  photon  travels  a  distance  p'  without  being 
absorbed,  defined  as 


d\>  e 
Avml 


rP ' 

.0  dD"  *;(r",t-p"/c) 


(14) 
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Specifically,  if  K^(r,t) 
relation  holds, 


K^(r',t')  then  the  following 


Gjnl  (  r '  / 1 '  ;  r ,  t ) 


-  1 


4n 


15) 


Physically,  the  requirement  is  that  over  the  mean  free  path 
of  an  emitted  photon  the  absorption  coefficient  does  not 
change  much,  i.e.  the  number  of  absorbers  encountered  during 
transit  remains  nearly  constant, 
as 


T  i ( p ' ; t )  can  be  written 


-  K'(r,t)  p' 

T.(p',t)  "  dv  <t>v  e  v  (16) 

'"’ml 

g 

Eq.  (16)  was  initially  introduced  by  Holstein  and  evaluated 
for  various  line  profiles. 

The  problem  is  now  reduced  to  the  following  steps:  (1) 
For  the  transitions  included  in  the  model  choose  an 
appropriate  line  shape  and  evaluate  Eq.  (16)  either  exactly 
or  approximately.  (2)  For  the  desired  geometry  obtain  the 
set  {G}  from  Eq.  (15).  (3)  Use  these  to  evaluate  the 

absorption  terms  in  Eq.  (12).  (4)  Solve  the  rate  equations. 

The  procedure  must  be  done  self-consistently  since  K  ' 
depends  on  the  solution  to  the  rate  equations. 


III.  MODEL 


A.  Nitrogen  Plasma 

Our  interest  is  in  studying  radiative  transfer  in  a  hot 

channel,  in  particular,  for  electron  temperatures  of  1.0  - 

3.0  eV.  In  this  temperature  range  the  constituents  are 

primarily  atomic  or  ionic.  We  have  included  in  the 

calculations  the  lowest  13  levels  of  N,  the  lowest  17  levels 

of  N+,  and  2  representative  levels  of  N++ .  Details  of  the 

13  14 

model  have  been  presented  elsewhere  '  .  All  bound-bound 

radiative  transitions  allowed  by  this  energy  scheme  are 
included  and  form  the  basis  for  the  line  radiation  discussed 
previously.  In  this  manner  the  dominant  uv  and  visible 
spectral  lines  are  transferred. 


B.  T(p) 

The  absorption  coefficients,  defined  in  Eqs.  (8)  and 
(9),  depend  on  the  line  shape,  We  have  assumed  this 

line  shape  is  given  by  a  Lorentzian,  i.e. 


*v,ml 


1  Av 

n  (v  -  vml)2  +  (Av)2 


(17) 


The  line  width  is  Av.  In  the  temperature  and  electron 
density  regimes  of  interest,  Ne  =  (  1016  -  1019  cm-3  ), 
Doppler  and  Stark  broadening  are  the  two  dominant 
mechanisms.  The  Lorentz  profile,  Eq .  (17),  is 
characteristic  of  the  Stark  broadening,  where  isolated 
spectral  lines  of  atoms  in  dense  plasmas  are  broadened  by 
electron  impact^.  However,  in  our  computations  we  consider 
the  effects  of  both  Doppler  and  Stark  broadening  whose 
actual  line  shape  should  be  a  Voight  profile.  In  the 
computations  described  below,  we  calculate  the  line 
broadening  of  both  as  a  function  of  frequency,  electron 
temperature,  and  electron  density  and  choose  Av  to  be  the 

greater  of  the  two.  The  Stark  effect  generally  dominates. 

1 8 

Additional  details  have  been  presented  elsewhere  . 
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For  the  Lorentzian  line  shape  Holstein  obtained  the 
asymptotic  functional  form  of  Eq.  (16)  as 


Tml(’>  - 
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ml 
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k  ^  is  the  absorption  coefficient  at  line  center.  The  exact 
integration  of  T(p)  is  shown  in  Fig.  2  along  with  two 
empirical  approximations.  For  x  ■  Kmlp  dotted  line 
corresponds  to 
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(20) 


Eq.  (20)  is  used  in  this  study. 

C.  Geometry  and  Gml 

The  results  discussed  in  section  IV  are  for  a  geometry 

where  radiation  is  assumed  to  flow  in  only  two  directions, 

but  along  one  dimension,  i.e.  r  =  r.  The  plasma  is 

contained  within  a  finite  region,  -L  <  r  <  L,  so  in  addition 

to  transfer,  energy  is  radiated  beyond  the  boundary.  For 

2 

this  geometry  the  4np'  factor  in  Eq.  (15)  is  excluded. 

When  evaluating  Eq.  (12)  care  must  be  taken  to  properly 
consider  three  possible  cases  which  depend  on  the  absorption 
length  at  line  center,  ,  and  the  distance  to  the  boundary 
edge,  L.  The  radiation  may  be  thick  in  both  directions, 

L-r  >  1/k  ^  and  L+r  >  1/k  It  may  be  thin  in  both 

directions,  L-r  <  1/Km^  and  L+r  <  1/k^ ,  or  thick  in  one 
direction  and  thin  in  the  other,  L+x  >  1/Km^  an(3 
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w 


L-x  >  1/Km^ .  Symmetry  considerations  enable  us  to  reduce 
computations  to  the  half  space  r  -  [ 0 ,  L ] .  A  typical 
absorption  term  for  the  case  where  the  line  is  thick  in  each 
direction  is 


6  J  0 


dp  N*_1(p,t) 


ml  Jl/K 


+  contribution  from  the  left 


Note  that  k 


<ml(c,t)  and  pm 


L-r.  Similiar  results  are 


obtained  for  the  other  cases. 


D.  Temperature  Equation 

As  radiation  is  transferred  from  one  local  region 
within  the  plasma  to  another,  cooling  or  heating  may  take 
place.  As  radiation  escapes  the  plasma  entirely,  a  net 
cooling  occurs.  To  model  this  phenomenon  the  relevant 
temperatures,  generally,  electron  and  gas,  must  be 
monitored.  For  hot  channels  created  by  intense  electron 
beams  the  plasma  electron  temperature  and  the  gas  kinetic 
temperature  frequently  reach  an  identical  value,  afterwhich 
the  subsequent  evolution  is  well  described  by  a  single 
temperature.  In  the  present  model  it  is  assumed  that  for  a 
given  set  of  initial  conditions  the  plasma  is  characterized 
by  a  single  temperature,  i.e.  T(r,t)  -  Te(r,t)  -  Tg(r,t). 
The  time  evolution  of  T(r,t)  is  obtained  from  the  energy 
conservation  equation, 

|  T(r,t)  (  Ne ( r , t )  +  Ng(r,t)  )  -  E  (r,0)  +  E  ( r  ,  0  )  + 


Ec(r,0)  -  Erad( r , t)  -  Ec(r,t)  . 


V-  \  W  V  A  - 


Here,  Ep(r,0),  E  (r,0),  and  Ec(r,0)  are  the  plasma,  gas 
kinetic,  and  chemical  energies  at  t  -  0  for  point  r, 
respectively.  Ne(r,t)  and  Ng(r,t)  are  the  electron  and 
total  gas  particle  densities  at  time  t  and  point  r.  E  (r,t) 
is  the  chemical  energy  at  time  t  and  Erad(r,t)  is  the  amount 
of  energy  radiated  from  (or  to  if  Erad  <  0)  point  r  at  time 
t.  Erad(r,t)  is  obtained  by  solving  the  following  equation, 

jr  Erad( r, t)  =  I 

m,  1 ,  z 

+  continuum  radiation  .  (23) 


V2-1 

ml 


hv 


ml 


Continuum  radiation  results  from  free-free  and  free-bound 
(radiative  recombination)  transitions.  In  this  model 
radiation  from  free-free  transitions  is  neglected,  uv 
radiation  from  all  2-body  recombination  events  is  assumed  to 
be  optically  thick  while  visible  radiation  from  these  events 
is  optically  thin.  Eqs.  (11),  (22),  and  (23)  are  the 
essence  of  the  model.  The  detailed  radiative  transfer  is 
contained  in  Eqs.  (11)  and  (23).  Eq.(22)  insures  energy 

conservation . 


IV.  NUMERICAL  RESULTS 


Numerical  integration  of  Eqs.  (11)  and  (23)  is 

accomplished  by  using  CHEMEQ ,  a  routine  designed  to  solve 
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stiff  ordinary  differential  equations  .  Symmetry  enables 

us  to  carry  out  the  integration  over  half  the  grid,  i.e. 

[0,L].  In  this  section  the  following  initial  conditions  are 

assumed:  The  size  of  the  region  is  1  cm,  i.e.  L  -  .5  cm. 

Initial  population  densities  are  distributed  uniformly 

across  the  grid  [ 0 , L ] .  The  values  are  N( 1 )+N( 2 )+N( 3 )  - 

3.715  x  1019  cm-3,  Ne  -  N+ ( 1 ) +N+ ( 2 ) +N+ ( 3 )  -  5.0  x  1016  cm'3, 

4  -3 

and  all  others  are  set  to  0.0  (  actually,  1  x  10  cm  ). 

The  specific  ground  N  and  N+  states  are  populated  according 
to  statistical  weights.  Two  initial  temperature  profiles 
are  assumed.  The  first  is  a  uniform  distribution;  Te  -  3.0 
eV  for  0.0  <  r  <  .5  cm.  The  second  is  non-uniform;  Te  -  3.0 
eV  for  0.0  <  r  <. .2  cm  and  decreases  linearly  to  a  value  of 
1.5  eV  at  r  -  .5  cm.  Setting  the  initial  population 
densities  and  temperatures  is  equivalent  to  setting  the 
initial  plasma  electron,  gas  kinetic,  and  chemical  energy 
distributions . 

Eqs.  (11)  and  (23)  are  integrated  until  an  equilibrium 
is  nearly  established.  During  the  integration  the  radiation 
terms,  for  example,  Eq.  (21),  are  updated  using  the 
solutions  to  the  rate  equations.  Initially  this  is  done 
every  10~10  sec  and  relaxed  to  1G~7  sec  at  longer  times. 
Figs.  3a  -  10a  correspond  to  the  uniform  temperature 
distribution  results,  while  Figs.  3b  -  10b  are  for  the  non- 
uniform  case. 

A.  Rate  Equation  Results 

Figs.  3  and  4  show  radiative  transfer  and  cooling  as 
evidenced  by  the  evolution  of  T(r,t)  up  to  t  *  10  3  sec  for 
r  -  0.0,  0.25,  and  0.5  cm.  For  both  cases  the  species 
rapidly  approach  a  local  equilibrium  while  the  energy  in  the 
system  readjusts.  This  occurs  because  the  initial 
conditions  correspond  to  a  nonequilibrium  situation.  The 
initial  relaxation  takes  approximately  10  ns  during  which 
time  radiative  transfer  has  not  begun  to  play  a  significant 
role.  For  the  uniform  case  cooling  is  nearly  constant 


across  the  grid  for  10  y sec.  The  outer  boundary  does  show  a 
lower  temperature  than  the  interior,  particularly  for  5-8 
fj sec.  This  is  due  to  more  energy  loss.  The  non-uniform  case 
shows  the  center  cooling  more  rapidly  than  the  boundary 
region.  Even  at  10  sec  a  temperature  gradient  remains. 

Fig.  5  shows  the  time  evolution  of  the  electron  density 
for  both  cases.  While  the  electron  density  was  initially 
constant  across  the  grid  for  the  uniform  case  the  extra 
energy  lost  from  the  boundary  lead  to  a  decrease  in 
temperature,  mentioned  above,  which  in  turn  lead  to  a  lower 
Ne  at  the  edge.  By  contrast  Ne  at  the  boundary  for  the  non- 
uniform  case  decreases  only  slightly  from  its  maximum  while 
the  interior  values  decrease  more  rapidly,  reflecting  the 
behavior  of  T( r , t ) . 

B.  Energy  Loss  Results 

The  time  histories  of  the  energy  profiles  for  r  -  0.0 
and  0.5  cm  are  presented  in  Figs.  6  and  7.  In  all  cases  the 
amount  of  chemical  energy  exceeds  the  gas  kinetic  energy 
which  is  greater  than  the  plasma  electron  energy.  Figs.  6a 
and  7a  show  the  extent  to  which  more  energy  is  radiated  from 
the  boundary  than  from  the  center.  Since  there  was  less 
energy  at  the  grid  edge  than  at  the  center  for  the  non- 
uniform  case,  the  opposite  situation  is  seen  in  Figs.  6b  and 
7b.  In  terms  of  energetics  these  figures  show  that  the  time 
scale  for  radiative  transfer  to  contribute  significant 
energy  loss  is  on  the  order  of  microseconds.  The  radiative 
energy  loss  values,  Erad(r,t),  for  r  -  0.0,  0.25,  and  0.5  cm 
are  show  in  Fig.  8.  By  10  //sec  these  terms  have  nearly 
saturated  as  have  the  energy  values. 

Throughout  the  integration  the  photon  energy  that  is 

created,  re-absorbed,  and  radiated  from  the  system  is 

monitored.  One  quantity  of  interest  is  the  radiated  energy 

-2  -1 

flux  for  each  line,  Qm^ »  given  in  units  of  eV-cra  -sec 
The  flux  radiating  from  one  side  is  calculated  according  to 


dr  Aml 1  hvml 


Table  1  shows  the  flux  values  for  the  major  contributing 
lines  as  a  function  of  time.  The  surge  in  emitted  radiation 
at  1  ns  corresponds  to  the  rapid  population  of  exited  N  and 
N+  states  during  the  initial  energy  redistribution  period. 
Subsequent  time  behavior  shows  a  slowing  down  of  the  cooling 
rate  as  radiation  energy  is  lost  from  the  system. 

Generally,  the  energy  loss  is  dominated  by  emissions  from 
the  nitrogen  uv  lines.  The  N+  uv  lines  and  N  optical  lines 
also  contribute  to  the  energy  loss,  though  to  a  lesser 
degree;  their  contribution  is  not  shown.  The  effect  of  the 
N+  visible  lines  is  negligible.  This  is  consistent  with  the 
corresponding  temperature  values. 


TABLE  1 

-2  -1 

Radiated  Energy  Flux  ( eV-cm  -sec  ) 


Time 

1493.3 

O 

Transition  Line  (A) 

1134.6  1199.9  1743.6 

1243.3 

Total 

100  ps 

7 . 8xl017 

1 . OxlO17 

2 . 2xl017 

3 . 9xl017 

2 . 3  xl 01 7 

1 . 9xl018 

1  ns 

8 .OxlO23 

6 . 3xl023 
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4 . 6xl023 

1 . 7xl023 

2.8xl024 

10  ns 

4. 4xl023 

3 . 6xl023 

3 . 4xl023 

2 . 6xl023 

2 . lxlO23 

2. OxlO24 

100  ns 

2  3 

3 . 9xl04J 

3 . 6xl023 

2 . 9xl023 

2 . 3xl023 

2. OxlO23 

1 . 9x1 02  4 

1  u  s 

3 . 4x10  2  3 

2 . 9xl023 

2 . 5xl023 

23 

2. OxlO4 

1 . 7  xl 0  2  3 

1 . 6xl024 

10  ps 

22 

9 . 2xl044 

6 . 7xl022 

22 

6 . 4xl044 

6 . 2xl022 

3 . 5x1 022 

3 . 9  xlO2  3 
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Line  Results 


It  is  convenient  to  define  a  line-integrated  mean 
z-1 

intensity,  Jm^  ,  as 


J^l(r't>  -J 

‘"'ml 


(25) 


In  terms  of  the  appropriate  absorption  term  and  population 
densities  it  follows  that 
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ml 
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ml 


N 


Z-1 
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.  (26) 


Pigs.  9-10  show  the  line  integrated  mean  intensities 

0  0 

for  2  transition  lines  (1199.9  A,  and  8211.8  A)  at  different 
times  (10  ns,  100  ns,  1/js,  and  10/us).  The  distribution  of 
radiation  for  these  lines  is,  therefore,  monitored  as  a 
function  of  time  and  space.  The  results  are  consistent  with 
the  associated  temperature  profiles.  As  expected  the  line 
intensities  for  the  uniform  case  show  much  less  variation 
over  the  grid  than  the  non-uniform  case,  particularly  at 
early  times.  As  radiation  transfer  occurs  energy  diffuses 
outward.  Fig.  9b,  in  particular,  shows  the  extent  to  which 
the  line  intensity  in  the  central  region  decreases  compared 
to  the  outer  region  as  a  function  of  time. 

D.  Discussion 

Diffusion  and  cooling  are  two  aspects  of  the  radiative 
transfer  process.  When  energy  is  distributed  uniformly 
across  a  confined  region  it  is  seen  that  the  outer  boundary 
cools  more  rapidly  than  the  interior,  thus  establishing  a 
temperature  (or  energy)  gradient  which  drives  subsequent 
transfer  and  net  cooling.  Radiative  transfer  initiated  by 
the  existence  of  a  temperature  gradient  behaves  differently. 


In  this  case  energy  diffuses  outward  from  the  outset  and  the 
hotter  inner  region  cools  more  rapidly  than  the  cooler  outer 
region.  Using  the  model  descibed  previously  the  above 
numerical  results  for  the  given  set  of  initial  conditions 
demonstrate  that  the  radiative  transfer  time  is  on  the  order 
of  microseconds.  This  can  be  understood  within  the  context 
of  the  following  simple  picture: 

If  radiative  transfer  is  considered  to  be  a  sequence  of 
photon  emissions  and  absorptions  it  should  be  possible  to 
estimate  bounds  on  the  time  required  for  radiation  to 
diffuse  from  the  center  of  a  region  to  a  boundary.  For  a 
specific  transition  line,  m*l ,  the  diffusion  time  may  be 
estimated  as 


t 


ml 


L  1 

(27) 

xml  Aml 

(28) 

where  L  is  the  distance  to  the  boundary  (.5  cm  in  this 
study)  and  is  the  "mean  free  path".  This  definition 

corresponds  to  the  following  picture:  After  photon  travels 


a  distance  it  is  absorbed  and  re-emitted  in  1/A 


ml 


the  effective  emission 
defined  as  A_,  -  Aml^Tml  w^ere  T-1  the 


ml 


ml 


(seconds).  It  continues  on  in  this  manner  until  reaching 
the  boundary.  In  Eq.  (28 
coefficient  is 
optical  depth. 

Within  the  context  of  our  radiative  transfer  model,  the 
definition  of  x^  is  not  obvious.  The  simplest  choice  is  to 
define  xm^  as  the  inverse  of  the  absorption  coefficient  at 
line  center,  x^  ■  1/k^.  T^is  puts  the  emphasis  for 
radiative  transfer  on  photons  that  are  emitted  directly  at 
the  center  of  the  frequency  profile.  The  probability  that  a 
photon  travels  this  far  without  being  absorbed,  as  given  by 
Eq .  (20),  is  0.67.  In  the  simple  diffusion  picture  it  is 
assumed  that  the  photon  is  then  absorbed.  This  gives  an 
upper  bound  on  the  diffusion  time  scale.  For  comparison, 
xm^  is  defined  as  the  distance  at  which  I’m]_  (  p  )  *  0.05.  In 
this  manner  absorption  is  guaranteed  after  the  photon  has 


18 


travelled  x  ^  -  177. 8/k^  .  The  later  definition  places  more 
emphasis  for  radiative  transfer  on  photons  emitted  from  the 
line  wings  since  these  are  the  photons  that  contribute  to 
the  asymptotic  behavior  of  T  ^ ( P )  •  This  gives  a  lower  bound 
for  the  onset  of  radiative  transfer. 

Table  2  shows  a  comparison  of  these  diffusion  times  for 
the  top  12  lines  contributing  to  radiation,  as  determined  by 


^ml 


at  t  «  10  ns.  The  diffusion  time  using  x 


ml 


1/K 


ml 


is 


i  w i  no 

denoted  as  tm1  and  tm-,  ^  is  the  estimate  for  the  larger 
value  of  x 


ml  “““  ml 

These  estimates  are  in  fact  consistent  with 


ml ' 

the  results  presented  in  Figs.  3  -  10. 


TABLE  2 


Diffusion  Time  Estimates  (sec) 
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CONCLUDING  REMARKS 


A  radiative  tranport  model  has  been  developed  which 
emphasizes  the  transfer  of  line  radiation  in  a  hot  nitrogen 
plasma.  The  radiative  transfer  equations  are  solved  self- 
consistently  with  the  corresponding  time-dependent  chemistry 
equations.  This  approach  allows  for  the  transfer  of  many 
lines  and  does  not  require  a  state  of  LTE  to  exist. 

Numerical  results  for  two  simple  test  cases  show  the  utility 
of  the  model  in  describing  radiative  transfer  effects.  A 
following  report  will  present  the  results  of  using  this 
model  to  account  for  radiation  transDort  in  the  realistic 


a)  Uniform  initial  conditions,  b)  Non-uniform 
initial  conditions. 


initial  conditions. 


ns  (top  solid  line),  100  ns  (dotted  line),  1  v s 
shed  line),  and  10  fjs  (bottom  solid  line).  a) 
form  initial  conditions,  b)  Non-uniform  initial 


REFERENCES 


1.  S.  Chandrasekhar,  Radiative  Transfer,  Dover,  New  York 
(  1960)  . 

2.  A.  L.  Crosbie  and  R.  L.  Dougherty,  J.  Quant.  Spectrosc. 
Radiat.  Transfer  2_5,  551  (  1981). 

3.  A.  S.  Eddington,  The  Internal  Constitution  of  the  Stars , 
Cambridge  Univ.  Press  (1926).  Also  Dover,  New  York 

(  1959)  . 

4.  J.  L.  Kulander,  J.  Quant.  Spectrosc.  Radiat.  Transfer  6, 
31  (1966). 

5.  B.  Davison  (with  J.  B.  Sykes),  Neutron  Transport  Theory, 

Clarendon  Press,  Oxford  (1957). 

6.  M.  A.  Heaslet  and  R.  F.  Warming,  J.  Quant.  Spectrosc. 
Radiat.  Transfer  6,  751  (1966). 

7.  D.  R.  Churchill,  S.  A.  Hagstrom,  and  R.  K.  Landshoff,  J. 
Quant.  Spectrosc.  Radiat.  Transfer  4^,  291  (1964  ). 

8.  L.  Baker,  CHARTB  Multigroup  Transport  Package",  SAND-79- 
0386,  March  1979. 

9.  T.  Holstein,  Phys.  Rev.  1_2,  1212  (1947);  8_3,  1159 
(1951). 

10.  L.  M.  Biberman,  Pure  Appl .  Chem.  1 3 ,  393  (1966). 

11.  J.  P.  Apruzese,  J.  Davis,  D.  Duston,  and  K.  G.  Whitney, 
J.  Quant.  Spectrosc.  Radiat.  Transfer  23,  479  (1980). 

12.  H.  W.  Drawin,  Z.  Physik  225 ,  470  (1969). 

13.  R.  D.  Taylor  and  A.  W.  Ali,  "Recombination  and 
Ionization  in  a  Nitrogen  Plasma",  NRL  Memorandum  Report 
5594 ,  (1985).  (AD-A  155  426) 

14.  Ronald  D.  Taylor  and  A.  W.  Ali,  J.  Quant.  Spectrosc. 
Radiat.  Transfer  35 ,  213  (1986). 

15.  R.  D.  Taylor  and  A.  W.  Ali,  "A  Collisional-Radiative 
Model  and  Saha  Decrements  for  a  Nonequilibrium  Oxygen 
Plasma",  NRL  Memorandum  Report  5792  (1986).  (AD-A168  082) 

16.  Ronald  D.  Taylor  and  A.  W.  Ali,  "Saha  Decrements  and 
Collisional-Radiative  Recombination  and  Ionization 
Coefficients  for  a  Nonequilibrium  Oxygen  Plasma",  J. 
Quant.  Spectrosc.  Radiat.  Transfer  v.  25, 

p.  213,  1986. 

17.  H.  Griem,  Plasma  Spectroscopy,  McGraw-Hill,  New  York 
(  1964)  . 


3] 


R.  D.  Taylor  and  A.  W.  Ali,  "Total  and  Partial 
Absorption  Coefficients  for  a  Nitrogen  Plasma",  NRL 
Memorandum  Report  5607  (1985).  (AD-A  157  534) 


T.  R.  Young,  " CHEMEQ  -  A  Subroutine  for  Solving  Stiff 
Ordinary  Differential  Equations",  NRL  Memorandum  Report 
4091  (1980).  (AD-A  083-545) 


DISTRIBUTION  LIST 

Chief  of  Naval  Operations 

Washington,  DC  20350 

ATTN:  Or.  C.  F.  Sharn  (0P09873) 

U.  S.  Army  Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground,  Maryland  21005 
ATTN:  Dr.  Donald  Eccleshall  ( DRX8R-BM) 

Or,  Anand  Prakash 

Office  of  Under  Secretary  of  Defense 
Research  and  Engineering 
Room  3E1034 
The  Pentagon 
Washington,  DC  20301 
ATTN:  Mr.  John  M.  Bachkosky 

Office  of  Naval  Research 
800  North  Ouincy  Street 
Arlington,  VA  22217 
ATTN:  Or.  C.  W.  Roberson 

Chief  of  Naval  Material 
Office  of  Naval  Technology 
MAT-0712,  Room  503 
300  North  Ouincy  Street 
Arlington,  VA  22217 
ATTN:  Or.  Eli  Zimet 

Commander 

Naval  Sea  Systems  Command 
PMS-405 

Washington,  DC  20362 
ATTN:  CAPT  R .  Topping 

CDR  W.  Bassett 

Air  Force  Office  of  Scientific  Research 
Physical  and  Geophysical  Sciences 
Bolling  Air  Froce  Base 
Washington,  DC  20332 
ATTN:  CAPT  Henry  l.  Pugh,  Jr. 

Department  of  Energy 
Washington,  DC  20545 

ATTN:  Or.  Terry  F.  Godlove  (ER20:GTN,  High  Energy  and  Nuclear  Physics) 

Dr.  James  E.  Lei.ss  (G-256) 

Mr.  Gerald  J.  Peters  (G-256) 

Joint  Institute  for  Laboratory  Astrophysics 
National  Bureau  of  Standards  and 
University  of  Colorado 
Boulder,  CO  80309 
ATTN:  Or.  Arthur  V.  Phelps 


33 


s  t  •"  «*  *  ,•  **  *  *■  *  n  ^ f.  ^  ^  ^  k  ^  k  ^  •  .k"'  k  *  I 


c-v-v- 


T 


y  ■  u  ■  T  ■  t 


r»  v*  v  *  4 


Lawrence  Berkeley  Laboratory 
jmversity  of  Califor'na 
Berkeley,  CA  94720 
ATTN:  Or.  Edward  P.  Lee 

Ballistic  Missile  Defense  Advanced  Technology  Center 

P.0.  Box  1500 

-*untsvi  lie,  AL  35807 

AT’N :  Or.  M .  Hawie  73MDS AT C- 1 ) 


Intelcom  Bad  Teen. 
p.O.  Box  8108 7 
San  Oiego,  CA  92138 
ATTN:  Or.  W.  Selpn 


Lawrence  Livermore  National  Laboratory 
University  of  California 
Livermore,  CA  94550 


ATTN: 


Or.  Richard  J.  Briggs 

Or.  Thomas  Fessenden 

Or.  Frank  Chambers 

Or.  James  W.-K.  Mark,  L-477 

Or.  william  Fawley 

Or .  Wi 1 1 i am  Bari etta 

Or.  Willi  am  Sharp 

Dr.  Daniel  S.  Prono 

Or.  John  K.  Boyd 

Or.  Kenneth  W.  Struve 

Dr.  John  Clark 

Or.  George  J.  Caporaso 

Dr.  William  E .  Marti n 

Dr.  Donald  Prosnitz 

Dr.  S.  YU 


Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93102 
ATTN:  Dr.  C.  Longmi re 

Or.  N.  Carron 


National  Bureau  of  Standards 
Gai thersburg,  MO  20760 
ATTN;  Or.  Mark  Wilson 


Science  Applications,  Inc. 
1200  Prospect  Street 
La  Jolla,  CA  92037 
ATTN:  Dr.  M.  P.  Fricke 

Or.  W.  A.  Woo  Ison 


I 

) 


Science  Applications,  Inc. 

5  Palo  Alto  Square,  Suite  200 
Palo  Alto,  CA  94304 
ATTN:  Or.  R.  R.  Johnston 

Or.  Leon  Feinstein 
Or.  Douglas  Keeley 

Science  Applications,  Inc. 
1651  Old  Meadow  Road 
McLean,  VA  22101 
ATTN:  Mr.  W.  Chadsey 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 


Spri  i 

ig. 

MD 

20910 

Mr. 

R. 

J. 

8iegal ski 

Dr. 

R. 

Cawl ey 

Dr. 

J. 

W. 

Forbes 

Dr. 

D. 

L. 

Love 

Dr. 

C. 

M. 

Huddleston 

Dr. 

0. 

E. 

Hudson 

Mr. 

W. 

M. 

Hi  nckl  ey 

Mr. 

N. 

E. 

Scofield 

Or. 

E. 

C. 

Whitman 

nr. 

M. 

H. 

Cha 

Or. 

H. 

S. 

Uhm 

Dr. 

R. 

Fiorito 

Or. 

H. 

C. 

Chen 

C.  S.  Draper  Laboratories 
Cambridge,  MA  02139 
ATTN:  Or.  E.  Olsson 

Dr.  L.  Matson 

Physcal  Dynamics,  Inc. 

P.0.  Box  1883 
La  Jolla,  CA  92038 
ATTN:  Dr.  K.  Brueckner 

Avco  Everett  Research  Laboratory 
2385  Revere  Beach  Pkwy 
Everett,  MA  02149 
ATTN:  Or.  R.  Patrick 

Dr.  Dennis  Reilly 

Dr.  D.  H.  Douglas-Hami lton 

Oefense  Technical  Information  Center 
Cameron  Station 
5010  Duke  Street 

Alexandria,  VA  22314  (2  copies) 


35 


Naval  Research  Laboratory 
Washington,  DC  20375 
ATTN:  M.  Lampe  •  Code  4792 

M.  Friedman  -  Code  4700.1 
J.  R.  Greig  -  Code  4763 

I.  M.  Vitkovitsky  -  Code  4701 

J.  B.  Aviles  -  Code  6650 
M.  Haftel  -  Code  6651 

T.  Coffey  -  Code  1001 

S.  Ossakow  -  Code  4700  (26  copies) 

P.  Sprangle  -  Code  4790 
Code  2628  -  20  copies 

A.  W.  All  -  Code  4700.1  (30  copies) 

D.  Book  -  Code  4040 

J.  Boris  -  Code  4040  > 

R.  Hubbard  -  Code  4790 

B.  Hui  -  Code  4790 

S.  Slinker  -  Code  4790 
G.  Joyce  -  Code  4790 
D.  Murphy  -  Code  4763 
A.  Robson  -  Code  4760 

D.  Colombant  -  Code  4790 
M.  Picone  -  Code  4040 
M.  Raleigh  -  Code  4763 
R.  Pechacek  -  Code  4763 
G.  Cooperstein  -  Code  4770 
Y.  Lau  -  Code  4790 
R.  Fernsler  -  Code  4790 
Code  1220  -  1  copy 
Records  -  1  copy 

Defense  Advanced  Research  Projects  Agency  * 

1400  Wilson  Blvd. 

Arlington,  VA  22209 
ATTN:  Dr.  S.  Shey 

Physics  International,  Inc. 

2700  Merced  Street 
San  Leanadro,  CA  94577 
ATTN:  Dr.  E.  Goldman 

Or.  H.  L.  Buchanon 
SDI0  -  DEW 

Office  of  Secretary  of  Defense 

Washington,  DC  20301 

ATTN:  Lt.  Col.  R.  L.  Gullickson 

Mission  Research  Corp. 

1720  Randolph  Road,  S.E.  . 

Albuquerque,  NM  87106 
j-'TTN:  Dr.  Brendan  Godfrey 

Dr.  Richard  Adler 
Dr.  Thomas  Hughes 
Dr.  Lawrence  Wright 

Director  of  Research 

U.S.  Naval  Academy 

Annapolis,  MD  21402  -  2  copies 


36 


Princeton  University 
Plasma  Physics  Laboratory 
Princeton,  NJ  08540 
ATTN:  Dr.  Francis  Perkins,  Jr. 

McDonnell  Douglas  Research  Laboratories 
Dept.  223,  Bldg.  33,  Level  45 
Box  516 

St.  Louis,  MO  63166 
ATTN:  Dr.  Evan  Rose 

Dr.  Carl  Leader 

Cornell  University 

Ithaca,  NY  14853 

ATTN:  Prof.  David  Hammer 

Sandia  National  Laboratory 
Albuquerque,  NM  87115 
ATTN:  Dr.  Bruce  Miller 

Or.  Barbara  Epstein 
Dr.  John  Freeman 
Dr.  John  E.  Brandenburg 
Or.  Gordon  T.  Leifeste 
Dr.  Carl  A.  Ekdahl,  Jr. 

Dr.  Gerald  N.  Hays 

Dr.  James  Chang 

Dr.  Michael  G.  Mazerakis 

University  of  California 
Physics  Department 
Irvine,  CA  92664 
ATTN:  Dr.  Gregory  Benford 

Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117 
ATTN:  0.  Straw  (AFWL/NTYP ) 

C.  Clark  (AFWL/NTYP) 

W.  Baker  (AFWL/NTYP) 

D.  Dietz  (AFWL/NTYP) 

Lt  Col  J.  Head 

Pulse  Sciences,  Inc. 

14796  Wicks  Blvd. 

San  Leandro,  CA  94577 
ATTN:  Dr.  Sidney  Putnam 

Dr.  John  Bayless 

Los  Alamos  National  Scientific  Laboratory 

P.0.  Box  1663 

Los  Alamos,  NM  87545 

ATTN:  Dr.  L.  Thode 

Dr.  A.  8.  Newberger,  X - 3 ,  MS-608 
Dr.  M.  A.  Mostrom,  MS-608 
Or.  T.  P.  Starke.  MS-942 


Institute  for  Fusion  Studies 
University  of  Texas  at  Austin 
RLM  11.218 
Austin,  TX  78712 

ATTN:  Prof.  Marshall  N.  Rosenbluth 

University  of  Michigan 
Dept,  of  Nuclear  Engineering 
Ann  Arbor,  MI  48109 
ATTN:  Prof.  Terry  Kammash 

Prof.  R.  Gilgenbach 

Directed  Technologies,  Inc. 

226  Potomac  School  Road 
McLean,  VA  22101 
ATTN:  Dr.  Ira  F.  Kuhn 

Dr.  Nancy  Chesser 

Titan  Systems,  Inc. 

8950  Villa  La  Jolla  Orive-Suite  2232 
La  Jolla,  CA  92037 
ATTN:  Dr.  R.  M.  Oowe 


Lockheed  Palo  Alto  Laboratory 
3251  Hanover  Street 
Bldg.  203,  Dept.  52-11 
Palo  Alto,  CA  94304 
ATTN:  Or.  John  Siambis 

University  of  Maryland 
Physics  Department 
College  Park,  MD  20742 
ATTn:  Or.  Y.  C.  Lee 

Or.  C.  Grebogi 

Science  Applications,  Inc. 

1710  Goodridge  Or. 

McLean,  VA  22102 
ATTN:  Dr.  A.  Orobot 

Dr.  K.  Papadopoulos 

Dr.  Jonn  p.  Jackson 
Kaman  Sciences 

1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80933 

DAS  I  AC  -  DETIR 
Kaman  Tenoo 

25600  Huntington  Avenue,  Suite  500 
Alexandria,  VA  22303 
ATTN:  Mr.  F.  Wimenitz 


'A-V.PV1 


